Towards a statistical theory of solid dry friction 
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Wearless dry friction of an elastic block of weight N, driven by an external force F over a 
rigid substrate, is investigated. The slider and substrate surfaces are both microscopically rough, 
interacting via a repulsive potential that depends on the local overlap. The model reproduces 
Amontons's laws which state that the friction force is proportional to the normal loading force N 
and independent of the nominal surface area. In this model, the dynamic friction force decays for 
large velocities and approaches a finite static friction for small velocities if the surface profiles are 
self-affine on small length scales. 

PACS numbers: 46.30.Pa; 64.60.Ht 
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I. INTRODUCTION 



The physics of solid dry friction is an old and fasci- 
nating field. Yet, many quite fundamental problems are 
still subject of debate. The basic phenomenological facts, 
though, have since long been known as the Coulomb- 
Amontons's laws of friction: (i) The frictional force is 
independent of the size of the surfaces in contact, (ii) fric- 
tion is proportional to the normal load, and (iii) kinetic 
friction is not (or not much) dependent on the,-yelocity 
and typically lower than the static friction forced. 

A simple explanation for these laws arises from Bow- 
den and Tabor's adhesion theory of friction, in which 
plastic deformation of the surfaces accounts for Ibc load 
dependence of real contact area and friction forceao. Fur- 
thermore, plastic deformation leads to a logarithmic time 
dependence of the static friction andJogarithmic veloc- 
ity dependence of the kinetic frictionSu. Although plastic 
flow is assumed to yield the main contribution to solid 
friction, other mechanisms may play a role as well. In 
particular, it was noticed long ago that elastical multi- 
stability and hysteresis also gives rise to frictionQ. 

Recently, steps have been taken towards the under- 
standing of wearless friction as a collective phenomenon, 
dominated by the competition of pinning forces emerging 
from rough sujj£aces and bulk elasticity, neglecting plastic 
dcformations&LJ. However, a quantitative understanding 
is still lacking. Attempts in this direction are inspired by 
studies of .the depinning transition of dxuop charge den- 
sity wavesB, interfaces in random-jnedia&EEl and of vortex 
lines in type-II superconductorstil, where the behaviour 
near the depinning threshold force turned out to be a 
non-equilibrium critical phenomenon described by new 
universal critical exponents and scaling laws. It is tempt- 
ing to assume that friction is a related phenomenon— In- 
deed, in a recent investigation by Cule and HwaEJ, a 
bead-spring model for friction has been considered which 
exhibits a depinning transition of the universality class of 



interface depinning. Bead or block chain models do how- 
ever not account for Amontons's laws. The aim of the 
present paper is to study a simple statistical model where 
friction solely arises from hysteretic elastic response, and 
to find whether it is nevertheless capable of reproducing 
these fundamental laws. The situation we consider is the 
weak pinning limit, where elastic multistability arises as 
a collective effect. The opposite limit of strong pinning, 
where multistability already emerges on the local scale of 
single traps, has been considered recently by Caroli and 
NozieresQ. 



II. THE MODEL 

To be specific, we consider an elastic body of weight N 
and linear size L, which is pulled over a rigid substrate, 
cf. Fig. [IJ. The two surface profiles, separated by a mean 
distance d, are parameterized by scalar height functions 
Z(x) and h(x), respectively, x denotes the 2-dimensional 
position vector in the reference plane parallel to the sur- 
face and the substrate. For simplicity, both surfaces are 
assumed to have the same statistical properties: They in- 
dependently obey Gaussian distributions with mean zero, 
characterized by a short-range pair correlation function 



(h(x)/i(x')) = ft 2 fc((x-x')/CT), 



(1) 



where h defines the width of the substrate surface, a 
is the typical lateral corrugation length, k(0) = 1 and 
fc(x) w for |x| > 1. Correspondingly, P = (Z 2 (0)) 
describes the width of the slider surface. Unless other- 
wise stated, we will however assume throughout the pa- 
per that both surface profiles obey the same distribution, 
hence I = h. 

Short range correlations characterize a macroscopically 
flat surface. Self-affine surfaces, on the other hand, are 
characterized, in Fourier space, by a height-height cor- 
relator {h k h k <) ~ 5(k + k')|k|- 2C - 2 , where C is the 
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roughness exponent. Fracture surfaces, for instance, typ- 
ically have C = 0.6. . .0.9t3. The power law behaviour 
is usually cut off below some wave number 1 /an, so it is 
reasonable to restrict the study first to short-range cor- 
relations. 



v(z) = v z n e(z), 



(5) 




FIG. 1. Cartoon of the model. A block of weight N is 
pulled over a rigid substrate. The two adjacent surface pro- 
files h and I, separated by a mean distance d, are enlarged. 
Typical overlap areas are shown. 

The lateral elastic properties of the bottom slider sur- 
face depend on the shape of the body and are usually 
quite complicated. In general, the elastic energy can be 
written asli-H 



d 2 x d 2 x'j a p 1 s(x - x!)d a rp(x)dL ra(x') 



(2) 



where r(x) denotes the local lateral displacement from 
the equilibrium position. For simplicity, we restrict our- 
selves to 7 Q( 3 7 i(x) = 7 (x) 5 ai 5ps and consider only two 
limiting cases: If the slider is a 2-dirpensional object, like 
a latex membrane pulled over a rodE-3, 7(x) = 7<5(x), i.e. 
the elastic interaction is local. If, on the other hand, 
the slider is a semi-infinite 3-dimensional object, elastic 
response is mediated by bulk elasticity and nonlocal in 
space with 7(x) 7/|x|. Both cases can be treated si- 
multaneously by introducing an exponent a, so that the 
elastic kernel in Fourier space scales like 



7(k) 



7 (x) ~ k a ~ 



(3) 



with a = 1 and 2 for bulk and surface elasticity, re- 
spectively. Consequently, the dispersion of the elastic 
energy will behave as k a . Generalizing Eq. (|2|) to a D- 
dimensional surface, the real space elastic kernel scales 
like 7 (x) - |x| 2 --°- Q . 

In the direction perpendicular to the reference plane, 
the slider surface would strictly also have to be treated 
as elastic, interacting with the substrate via a hard wall 
potential (in the absence of adhesion forces). To make 
the model analytically amenable, however, we allow the 
surfaces to overlap and introduce a repulsive potential 
V(z) that depends on the local overlap 



z(x) = /i(x + r(x, t)) + Z(x) - d. 



(4) 



This potential is used to mimic vertical elasticity. We 
choose 



where Q(z) is the Heaviside step function, and n > 1. 
With n = 3/2, the Hertzian result on the distance depep. 
dence of the repulsive force between two elastic spherestil 
is reproduced (this nontrivial dependence being the result 
of the interplay between Hooke's law and the spherical 
geometry). As it turns out, the results do not depend 
very sensitively on the chosen value of n. Note that set- 
ting V(z) = for z < via the Heaviside function is 
natural and justified, in the absence of adhesion. 

Finally, the total driving force F is applied homoge- 
neously as a force density f = F/L 2 . This choice appears 
to be natural in the case of bulk elasticity, whilst for a 
membrane-like slider the external force should rather be 
exerted at one border. 

Assume that, with the slider lying at rest on the sub- 
strate, we turn on the external driving force F. If F is 
large enough, the slider will start to move, and the inter- 
action between the slider and the substrate will generate 
a dynamic friction force Fjv- = — F (as we will show be- 
low) which leads, after a transient time of acceleration, 
to a constant average velocity v of the center of mass of 
the slider. In the steady state, the equation of motion of 
a point r(x, t) of the slider surface can be written as 

r) (r(x, t) - v) = J dV 7 (x - x')V 2 r(x', t) + f (6) 



d 



dr(x,t) 



V[h(x + r(x,t))+l(x)-d]. 



The term on the left hand side accounts for surface 
phononic damping within the slider, with a damping co- 
efficient n. Since Eq. (^) is written in the laboratory 
frame, the center of mass velocity v has to be subtracted 
from the local velocity f . In this way we make sure that 
rigid sliding, i.e. r(x) = v£ with a constant velocity v, 
is frictionless. The neglect of an inertial term is jus- 
tified by our primary interest in the depinning region, 
where kinetic energy is relatively small. If we started, on 
the other hand, with a theory including an inertial term 
only, viscous friction would nevertheless be generated by 
the nonlinear random force after eliminating short wave- 
length displacement modes. 

The mean separation d between slider and substrate 
obeys a similar equation of motion 



-N - 



d 2 x ^V[h(x + r(x, t)) + i(x) - d], 
od 

(7) 



which has to be solved simultaneously with Eq. (|g). 
Again, overdamped motion is assumed with another fric- 
tion constant. In the analysis we will take A — > 00, 
because in the thermodynamic limit of infinitely large 
system size L, fluctuations of the mean distance d will 
vanish, so that it can be treated as a constant parameter 
that has to be determined self-consistently. 
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Simple scaling analysis shows that the model is domi- 
nated by two dimensionless quantities: 



£ = Vo^/ja 2 



and 



TV = N/V h n 



(8) 



where £ is the ratio between typical overlap and elastic 
forces on short scales. In this paper, we will restrict our- 
selves to £ <C 1, corresponding to weak disorder or weak 
pinning, and 1 <C Af <C (L/cr) 2 . Under these conditions, 
Af will turn out to be of the order of the number of con- 
tact points between the two surfaces (Vb/i n_1 (7 2 is the 
typical overlap force at a single contact), hence the lat- 
ter condition ensures that the real contact area is smaller 
than the nominal area L 2 and bigger than the typical size 
a 2 of one point of contact. Note that Af <C (L/a) 2 can 
be rewritten as p <C po , where p = N/ L 2 is the nominal 
pressure and po — Voft™ -1 a typical pressure at contact 
points. 



is weak, influencing only the strength of the logarithmic 
correction. 

The proportionality between load and contact area is a 
generic feature in our model, as long as the height prob- 
ability distribution P(h) decreases at least exponentially 
for heights \h\ 3> h, where h is the surface profile width. 
The only restriction imposed on the choice of the poten- 
tial V(z) is that it is cut off for z < 0, which should be 
true for any effective, adhesionless interaction, and that 
it increases slower than exponentially, or, in the case of 
a Gaussian height probability distribution, at most ex- 
ponentially; the exact form for z > does not matter. 
Under these assumptions, the expectation value (consid- 
ered to be a function of d) satisfies 

(V{h -d)) = / V(h - d)P(h)dh ~ P(d) (13) 



III. STATIC PROPERTIES 

Let us first consider the limit of vanishing surface elas- 
ticity £ — > (which implies that there is no friction force 
and hence no pinning at all) and zero external force F. 
In this limit r(x, t) = 0, and in Eq. (Q) the height profiles 
h and I can be averaged over. We find 



1 



V Q h n - x <T 2 



d 2 x (V'[h(x) + l(x)-d]) 



(9) 



^ > ■ -d 2 /4h 2 



<■„ I - J (2h/d) n e 



[i + o(Ch/d) 2 )), (io) 



where c„ = T(n + 1)/V4vr. This implicitly determines 
d as a function of Af and the surface roughness: d rs 
2/iln 1 / 2 (c n L 2 1 'No 2 ') . The real contact area 



A r = L 2 (<d(h{x) + l{x)-d)) 



(11) 



can be calculated in the same limit. Note that, while lat- 
eral elasticity is suppressed in the limit £ — > 0, the "soft" 
interaction potential V(z) that mimics vertical elastic- 
ity still allows for a finite contact area. Making use of 
Eq. (B), it is given by 



Ma 2 , 



47T C„ 



In" [c n L 2 /Na 2 ] 



(12) 



Apart from logarithmic corrections, A r is thus propor- 
tional to the normal load, independent of the total sur- 
face area 1 . The dependence on the potential parameter n 



to leading order. With the definition in Eq. (^Tj) and 
N/L 2 = (V'(h + I — d)), the proportionality between A r 
and N is a consequence of this result. 

We next assume that the elasticity constant £ is finite 
and set up a perturbation expansion in the local displace- 
ment field r(x, t) about an undistorted configuration, first 
with F — 0. For our further calculation it is useful to 
generalize the model to a Z?-dimensional slider surface, 
redefining Af correspondingly to Af = N/Voh n ~ 1 <j D . The 
definition of £ is unaffected. The set of equations of mo- 
tion (j^) and ([?]) can then be solved by iteration, order 
by order in the strength of the interaction potential Vo, 
using a diagrammatic technique introduced inlill. At each 
step, the random profiles h and I are averaged over. In 
the following analysis, the mean distance is assumed to 
be determined self-consistently by Eq. (^) in the limit 
A — > oo, i.e. it is treated as a constant. 

In order to set up the perturbation theory, we trans- 
form the equation of motion (|^) to Fourier space: 



Go '(k.w) ? k , w = (27r) D+1 6 D (k)S(u;)f 



(14) 



d D xdte-^-' lut 



d 



dr(x,t) 



V[h(x + r(x,t))+l(x)-d], 



where the bare propagator is given by 

G (k,^) = ( 7 |kr+i w )- 1 . (15) 

Expanding the last term in Eq. ( |l4| ) in powers of r (x, t) , 
one gets for the displacement field to first order in Vq 



, the proportionality between load and contact area was 
shown for a purely elastic surface with Gaussian distributed 
profile. In the model considered there, surface asperities were 
assumed to be spherically shaped and to obey the Hertzian 
distance-force relation. 
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fk,w = -Gro(k,a;) / d D xdt 



— ikx— iwt 



(16) 



v/i(x) y'[/i(x) + z(x) - rf] + o(y 2 ). 

For the displacement correlation function in lowest non- 
vanishing order follows 



(iWkw) = -(27r) D+2 d(Lo)S(Lu')d D (k + k') X 
|Go(k, w )| 2 



(17) 



d D xe-' lk * x 



V r| r=0 + ^ + *(*) - ^ WO) + - rf D ' 

ee C 2 (x + r,x;d) 

Since the surface height profiles are Gaussian distributed, 
the pair correlator C 2 will be a function of the first (which 
we have chosen to be zero) and second moments of the 
distribution only. More specifically, only the sum of the 
second moments enters the definition of C2, so we can 
define 

C 2 {K h (x 1 ) + #j(x 2 ); d) = C 2 (xi,x a ; d), (18) 

with K h (x) = (h(x)h(0)) and JQ(x) = (Z(x)Z(O)), where 
the average is taken over the respective height distribu- 
tions. The exact form of C% (K, d) is analyzed in ap- 
pendix |A|. The main result of this analysis is that the 
value of C2 is proportional to the normal load N if 
l x i|, l x 2| < c. For the real space displacement correla- 
tion function follows (with a <C |x| <C L) 



\{ [r(x) - r(0)] 2 > = Jjl - e ikx )(f k ^ f- k ,o> 



(19) 



1 - e 1 



kx 



k.p 



7 2| k |2a 



^i^x 2 p 2 e- ikX2+ip(x2 - Xl) C 2 (x 1 ,x 2 ;d) 



a 2 ME 2 (a/L) D x 



D c -D \ a 



. . D € -D 

HI for D < D c 



for D = D r 



where f k is short hand fo r f -£ ^yj - In the last equation, 
we have made use of Eq. ( |All| ) from App. |A|, with k » 1. 
The (upper) critical dimension is given by 



D c = 2a 



(20) 



implying that the physically relevant situation D = 2 
and a = 1 corresponds to the marginal case. Ignoring 
for a moment the dependence on length scales, the result 
is proportional to ME 2 ~ iVVb/7 2 ; the linear dependence 
on N stems from replacing Voe~ d ^ 4h (being the lead- 
ing dependence on the mean distance d) by N/L D via 
Eq.(§). 

Like in other cases of collective pinning, one can 
now determine a Larkin length Ll from the condition 
( [r(L L )-r(0)] 2 ) « er 2 E!l, i.e. L L is the length scale where 
typical displacements become of the order of the corru- 
gation length of the substrate. This yields 



L L ^a(L D e/a D M£ 2 ) 1 ^ 

= a[ 1 2 o 2( - 2 -^/V Q h n+1 p] 1 l\ 



(21) 



where e — D c — D. Note that Ll depends, apart from a 
combination of intrinsic model parameters, only on the 
nominal pressure p = N/L D . 

What are the consequences for the static friction force, 
which in our model is identified with the critical force at 
the depinning transitiont2l? Ll will be smaller than the 
system size only if we choose our model parameters such 
that Af ^> £~ 2 (L/a) D ~ £ . In this case, it is possibleio es- 
timate the static friction force in the standard wayE3iEl: 
On small length scales, the elastic energy is dominant 
compared to the interaction with the random potential, 
so that adjacent sites move coherently. On larger length 
scales, the random forces become more relevant until on 
the scale of the Larkin length Ll, elastic and pinning 
forces are of the same order of magnitude, so the elas- 
tic manifold becomes able to explore the inhomogeneous 
force field emerging from the randomly distributed con- 
tact points on scales larger than Ll- Regions of linear 
size Ll can hence be assumed to adapt independently 
to the disorder, each of them giving an independent con- 
tribution to the pinning force of the order of the elastic 
force on this scale. The force density in a Larkin region 
is thus given by /^(Ll) ~ jaL^ a , which is the product 
of a typical gradient a / X 2 , a typical value of the elastic 
kernel 7L^ £l+2 ~ a , and the area L® of the region. The 
total friction force is thus of the order 



Ffr « [LILlT Lt>ffr(L L ) « jL u aLl a , 



(22) 



which is proportional to the nominal area L D and hence 
violates Amontons's first law. 

With J\f(a/L) D < 1 and £ < 1, however, L L will typ- 
ically be much larger than L. The fact that in solid dry 
friction, the Larkin length is typically 3> -LrJias been re- 
cently remarkedjalso by Persson and TossatiQ and by Car- 
oli and NoziereaJ. In the following, we will restrict our- 
selves to this situation. In the physical situation D = 2 
with bulk elasticity a = 1, being just at the marginal 
dimension D c (i.e. e = 0), Ll is exponentially largelljCI: 



cexp c 



L L 



a D M£ 2 



(23) 



where c is a dimensionless constant. 

The elastic response of the surface also allows for a re- 
duction of the potential energy of the slider by decreas- 
ing the mean distance d. Technically, this appears in the 
form of a correction factor to the r.h.s. of Eq. (0). It is 
calculated by expanding V'[h(x+r(x)) +l(x) — d] in r(x) 
and substituting for r(x) the first order expression (p^), 
yielding 
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(V'[h(x + r(x)) + Z(x) - d]) = (V'[h(Q) + l{0) - d]) 



G o (k,0) J d D x'e-' l]sx ' 



V 2 r \ r=Q (V[h(x! + r) + Z(x') - d]V'[h(0) + 1(0) d]) 

+ O(V h n - 1 £ 2 ). (24) 

In this order of perturbation theory, the correction can 
be represented by a factor (1 — c£(h/d) n ) with another 
dimensionless constant c. Since this correction factor is 
smaller than 1, the self-consistent determination of the 
mean distance d via Eq. (^) will in turn lead to a slightly 
decreased value of d — ► d— Ad, so that the enlarged value 
of the leading factor 



-(d-Ad) 2 /4h 2 



(25) 



just compensates the reduction. A similar correction in 
Eq. (|l2]), together with the modified d, leads to a modi- 
fied contact area A r — + A r + SA r . 

In order to determine the sign of the relative correction 
8A r /A r , one has to compare the magnitude of the rela- 
tive corrections to (V'(h + I — d)) and to A r , as d efined 
in Eq. (|ll|), respectively. Using the result of Eq. (A14) 
in App. ^ one finds that the relative correction to A r is 
smaller by an approximate factor of T(n)T(n + 1)/T(2n) 
(which is smaller than 1 for n > 1). Consequently, the 
increase in A r due to — Ad < in the factor Eq. (E5T) will 
dominate, so the ratio 6A r /A r (which is of order £ ) is 
positive. 

If there is a friction force at all, it should depend on 5 A 
and vanish for 6 A = 0, because a completely rigid surface 
is never pinned. We can hence give a dimensional argu- 
ment for a characteristic friction force Ffr. The simplest 
way to estimate Fp. is to write it as the product of the 
excess contact area 6 A 
force density Voh n /a = 



a D Af£ and a typical lateral 

l-a g 



F, 



Jr 



7(7 



-D+1- 



a N£ 2 



(26) 



Here, a logarithmic correction of the order \n{L D / No~ D ) 
has been omitted. Defining as usual the friction coeffi- 
cient fi = Ffr/N, this expression corresponds to a value 
of \i of order £h/a which depends on the ratio between 
elastic and repulsive forces on small length scales but not 
on the load. Below, we will show that this estimate gives 
indeed the right order of magnitude of Ffr . 



IV. KINETIC AND STATIC FRICTION 

Next, we consider the case of a moving slider, v = |v| ^ 
0, driven by a finite force F. The perturbative expansion 
is now set up in 



u(x, t) = r(x, t) — vt, 



(27) 



relative to steady, rigid sliding. For further simplifi- 
cation, we allow only for displacements in the direc- 
tion en = F/|F|, corresponding to a further restriction 



7a/37(?(x) = 7(x)5 Q , 7 5 / 35<5 Q ,i in Eq. (g), where 1 denotes 
the direction en. Hence, u(x,i) = u(x, t)ey. It has been 
argued in a closely related context (the depinning of a 
driven flux line in a random medium) that this restric- 
tion will not alter ±he critical dynamics in the average 
direction of motionE3. This approximation will however 
overestimate the force needed to overcome a repulsive 
trap, because asperities cannot avoid each other by sim- 
ply bending away. Recently it has been shown explicitly 
that this avoiding process completely rules out local mul- 
tistability in the case of isolated isotropic traps, while 
multistability is re-established for anisotropic trapsE3. 
Full account of D-dimensional lateral elasticity has been 
taken in a simpler version of our model, where the ran- 
dom potential energy is proportional to the product|-e|f 
two charge densities on the two elastic manifolds, incJ; 
there, however, the case of Ll <C L was considered. 

To find the average dynamic friction force to lowest 
order in perturbation theory, we follow a procedure used 
by Feigel'man for driven interfacesEl He calculated the 
first correction to the mobility constant rj, treating the 
average velocity v as a parameter which has to be deter- 
mined self-consistently. Technically, we proceed as fol- 
lows: First, replace r(x,t) by vt + u(x,i) in Eq. (Ill), 
and expand the random potential part in this equation 
in powers of u(x, t). Now, insert the Fourier transform 
of Eq. (|l6|), having replaced the l.h.s. of this equation by 

u and the argument of h on the r.h.s. by x + vt, into 
the next order expansion term in Eq. (fijl), 



d x dt c 



— ikx— iut 



d 2 



(28) 



and perform the average over h and /. One finds to sec- 
ond order in Vq 



/= ^::' 2 " r(2 r" / ^.n,* 



d 2n+l 



(29) 



p,Sl 



d D xdte ipx+i 



Qi 



2h 2 



exp 



2h 2 + K h {yi + vt)+Ki{yi) 



Here, the representation of the random potential correla- 
tion function obtained in App. [a| with the dimensionless 
function /2(k) as defined in Eq. ( |A14 ), has been used. 

The r.h.s. of Eq. (p9[), multiplied with the full area 
Ld of the slider, is identified with the total friction force 
Ffriv) = L D ffi.(v). The resulting expression for Ffr(v) 
can be written in the form 



F^V) 



(30) 



C2(o-f,Cr(X- D) 



Lx \X 



2 ° + (7t) 2 £ 



2D 
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where C2(p, p'), defined in Eq. (A18) in App. [A|, is the 
Fourier transform of the random potential correlator de- 
fined in Eq. (17). X = a P an( i £ = crq are dimension- 
less integration variables, and we introduced the velocity 
scale 



V = 7CT 1 a /i], 



(31) 



which is a typical relaxation velocity. 

Written in this way, the integral on the r.h.s. of Eq. (|3C 
is just a number, the velocity dependence entering only 
via the ratio v/vq. Hence, we can write the friction force 
in the form 



(32) 



Ffr(v) ~Ffr4>[- 



where <fi(x) is a dimensionless function that depends on 
the explicit form of V(z) and the statistics of the surface 
profiles. A detailed discussion of the velocity dependence 
in 4>(x) is given in App. |^. Here, we will only summarize 
the results. 

At v w Vq, to begin with, Ffr(v) is indeed found to be 
of the order Ffr, which followed from our naive estimate 
Eq. (|2li|). The amplitude Ffr, implying in particular the 
proportionality between Ffr(v) and TV, follows for similar 
reasons as those that led to Eq. (|l9|), making use of the 
random potential correlator calculated in App. |A|. Thus, 
our dimensional argument ( |26| ) is justified a posteriori. 

In the high velocity regime d > »o, the friction force 
decays proportional to 1/u, independent of D and a. 
Note that as we have not included inertial terms in the 
equations of motion (|^) and (0) , the experimentally ob- 
served velocity strengthening behaviour of the friction 
force for high velocities can not be reproduced by our 
model. 

In the most interesting regime v <C vq, the behaviour 
crucially depends on the characterics of the surface profile 
correlator k(x) and on the dimensionality. If the height- 
height correlator k(x) is analytic in the origin, the friction 
force reaches its maximum (of order Ff r ) at v ~ and 
it decreases for v <C v as (v/vo) 1 " e ^ a - This is the usual 
contribution to friction from phononic damping, vanish- 
ing for v — ► if e/a < 1. In the marginal dimension 
D = 2, with bulk elasticity, Ffr(v) thus depends linearly 
on v (times a logarithm of v/vq, cf. Eq. (B4) in App. |b|) 
in the small velocity regime. 

In order to find a finite static friction force, however, 
Ffr(v) has to bend towards a non-zero value for v — > 0. In 
related problems like the depinning of driven interfaces 
and charge density waves, it has been shown that a finite 
pinning threshold appears due to contributions to 4>{x) on 
length scales larger than the Larkin lengthEl On these 
length scales, configurational multistability emerges as 
a collective effect, leading to collective pinning. This is 
reflected, in a renormalization group treatment, by the 
renormalized random force correlator developing a cusp- 
like singularity at the origin. Since in our case Ll is 



typically larger than the system size, the collective pin- 
ning mechanism is absent. We thus have to find criteria 
for the existence of local multistability. 

As mentioned above, real world surfaces often have 
self-affine properties spanning several orders of magni- 
tude. A (one dimensional) surface profile, for instance, 
that is short range correlated on the scale a, and self- 
affine with a roughness exponent £ = 1/2 below this 
scale, is described, in real space, by the correlator 



K h (x - x 1 ) = (h(x)h(x')) = h 2 e-\ x - x '\' a . 



(33) 



If the self-affinity covers the length scales from a down to 
a microscopic length scale a, (h(x)h(x')) will exhibit the 
cusp only when considering it on a coarser scale than a; 
on finer scales, it is analytic again. In Fourier space, the 
height-height correlator corresponding to such a profile 
is given by 



K h {q) 



2 e -(9 Q ) 2 
1 + M 2 



(34) 



A surface with this correlation function can be generated 
by an Ornstein-Uhlenbeck procesafJ: Given a stochas- 
tic spatial noise ((x) with (£(x)) = and (£(xX(x')) = 
(2h 2 / a)S a (x — x'), where 5 a produces short range corre- 
lations over a length scale a, the profile h(x) obeying the 
differential equation 



dh 1 , / \ »r \ 

= h x +C (&) 

ax a 



(35) 



has the desired statistical properties. 

For such surface profiles, we find that 4>(v/vo) takes a 
finite value of order 1 for voia/a) 1 ^ <C v -C v . In this 
case, the total friction force is indeed of the estimated 
order Ffr and almost constant for v < Vq, and decays 
~ 1/v for v ^> vq. More generally, a finite static fric- 
tion force is found if the (unrenormalized) surface profile 
correlator fc(x) has a cusp in the origin (precisely, we 
need lim a;|| _ > o + d x ..k{x) ^ 0, where xy is the component 
of x parallel to f), which corresponds to the first spa- 
tial derivative of k(x) undergoing a jump (of order 1) at 
x = 0. A surface characterized by such a correlator has 
local slopes that may take arbitrarily high values, even- 
tually leading to multistability even for arbitrarily small 
roughness h. 

Note that for a <C a, the lower velocity scale 
vo(a/a) 1 ^ a will typically be so small that it has no sig- 
nificance for a finite sample, since velocity fluctuations 
will become so large that the slider gets pinned. Conse- 
quently, for L finite and hence in any experimental situ- 
ation the regime v <C VQ(a/a) 1 ^ a where <j)(v/vo) <C 1 is 
unlikely to be observable. 



G 



V. STRONG PINNING 



where 1 <C N <C {L/a) 2 has been used. 



Our results have been obtained using the framework of 
perturbation theory about weak disorder, which is suffi- 
cient in the weak pinning limit and which is the nat- 
ural starting point for a renormalization group analy- 
sis that reveals configurational multistability on length 
scales larger than L^. In principle however, it is possi- 
ble to have multistability already locally on the scale of a 
single trap*-. This situation corresponds to the strong pin- 
ning limitEa, which cannot be treated successfully within 
finite order perturbation theory. 

Such a situation was recently considered by Caroli and 
Nozieres (CN)0. They consider two flat surfaces with 
a sparse distribution of bumps and sinks, where 'active' 
traps are formed when two adjacent asperities are in con- 
tact. Writing the interaction at an active trap as a poten- 
tial energy V(p) that depends on the distance p between 
their centers, they derive a criterion for the existence of 
local (single-site) multistability which reads 



VI. CONCLUSION 



m&x\d 2 p V(p)\ > Ea, 



(36) 



E is Young's modulus - which, in D — 2 and with a = I, 
is our 7 - and a a typical length scale of the trap. Similar 
considerations have been applied for instance in mean- 
field like descriptions of driven interfaces or charge den- 
sity wavesEH. The frictional force is then proportional 
to the typical energy gap at a spinodal jump, multiplied 
with the density of active traps. 

To make contact with the CN model, we can identify 
the overlapping asperities in our model (which are rela- 
tively few in the limit of small normal load) with these 
active traps. We aim to derive an estimate for the onset 
of local multistability in the case of short range corre- 
lated surface profiles, with a typical corrugation length 
<7, and a potential V(z) = Voz 2 6(z), in the physical sit- 
uation D = 2 with bulk elasticity (a = 1). A typical 
energy that can be stored in one active trap can roughly 
be estimated as Voh 2 a 2 , and the second derivative with 
respect to a lateral displacement will be of order Voh 2 . 
Comparing this with the elastic energy term 7c, we find 
as condition for the presence of local multistability that 



8 >0(1), 



(37) 



a condition which violates our assumptions. 

This restriction can be relaxed to some extent if one 
looks for the occurance of the first appearance of a mul- 
tistable site in a finite sample when tuning, for example, 
the potential amplitude Vq. The largest asperity in a 
sample of linear size L will have a height h nmx of order 
2h In 1 ' ' 2 ((L/a) 2 ), and with the expression for the mean 
distance d given after Eq. ([)]) follows h mayi — d « In W. 
Using further that the second derivative of a Gaussian 
correlated surface profile will typically take a value of or- 
der /i m ax/c 2 at this maximum, we arrive at the criterion 



£ lnA/"> 0(1) 



(38) 



Expressions (|26|) and ( J32| ) are the main results of this 
paper. They describe a friction force which depends lin- 
early on the weight N (up to logarithmic corrections), 
but not on the nominal contact area L D and hence ful- 
fills Amontons's laws. The dependence on the potential 
parameter n is weak. It is noteworthy that this result 
was obtained in the limit A r <C L D , where statistics are 
dominated by rare events. 

Trying other forms for the overlap interaction poten- 
tial V(z), it turned out that the crucial ingredient lead- 
ing to Ffi, ~ N is the cutoff below 2 = 0. If one aban- 
dons this cutoff and uses, e.g., an exponential potential 
V(z) = Vq e z / z ° , the proportionality is no longer valid; for 
the exponential potential, for example, one finds instead 
to leading order that Ffr ~ N 2 . 

To summarize, we have introduced a new stochastic 
model that incorporates the interplay between bulk elas- 
ticity and surface roughness in solid dry friction. To our 
knowledge, it is the first purely elastic model that treats 
solid dry friction as a collective phenomenon and repro- 
duces the correct load dependence of the friction force 
Ff r , known as the Coulomb- Amontons's laws. For high 
velocities v, Ff r decays like 1/v, while the behaviour for 
small v depends on the surface profile statistics: For 
a smooth surface, the static friction force vanishes in 
the weak pinning case, while it is finite if the surface 
is characterized by a non-analytical height-height corre- 
lator. We have also given an estimate for the appearence 
of local multistability in the case of smooth interfaces. 

We have not considered thermal effects in our study, 
because in most situations where temperature changes do 
not strongly affect the mechanical strengths of the sliding 
bodies, the friction coefficient is fouad to be basically in- 
sensitive to temperature variationsc3. An example where 
friction does however strongly depend on the temper- 
ature is the case of rubber sliding over hard surfaces, 
where the friction properties are intimately connected 
with the temperature dependent visco-elastic properties 
of the rubber material^. In that specific case, the fric- 
tion coefficient r\ in our model would become temperature 
dependent, modifying the velocity scale (cf. Eq. (31)) 
while leaving the friction force amplitude Ff r (Eq. (26)) 
constant. 

The present study can be extended in many directions. 
First, let us shortly focus on the damping term in Eq. (jg). 
In general, also ry is non-local in space (and time): For 
long wavelengths, it behaves like r}(k) ~ k s , where the 
precise value of s depends on the damping mechanism 
under consideration. Attenuation rates of surface waves 
for a semi-infinite elastic body with a rough surface lead 
to s = 3 or 4 in the long wavelength regimcEa. In this 
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case, the dynamics are dominated by a modified dynam- 
ical critical dimension D c d = 2a — s, while the static 
properties are unaffected. 

In a situation where the Larkin length is smaller than 
the system size, the perturbative results hold only in the 
limit of large velocities; for this situation, the critical 
dynamics close to the threshold force remain to be ana- 
lyzed. Finally, inertial terms can be included to properly 
describe the high velocity regime. 

It is a pleasure to thank Jan Kierfeld, Tim Newman 
and Kirill Samokhin for useful discussions. This work has 
been supported by the German Israeli Foundation (GIF). 



APPENDIX A: RANDOM POTENTIAL 
CORRELATION FUNCTION 



Let v n (z) = Q(z)z n , cf. Eq. 
the pair correlator 



We want to calculate 



C 2 (K h { X ) + Ki{ X ');d)= (Al) 
(V[fc(x) + J(x') - d]V[h(0) + 1(0) - d]), 

cf. Eq. (!?]), for fixed x and x'. The first step is to rewrite 
Eq. ( Al) using the Fourier transform of V(z) = J q e lqz V q , 
followed by performing the disorder average over the 
Gaussian height fields h and I. This results in 

q,q' 

(A2) 

V q and Vgi are Fourier transformed back to V(z) and 
V(z'), and the Gaussian integral over q and q' is carried 
out. Before writing it down, we go over to dimensionless 
functions and variables. We choose 



l = y/(h?) + (P) 



(A3) 



as the length scale in the direction perpendicular to the 
surface plane. Let 



d _ K h (x) + K^) 

1 ' K — 12 ■ 



(A4) 



so k ranges between and 1, with k = 1 for x = x' = 
and k — > in the opposite limit |x|, |x'| ^> a. Finally, let 
us define the dimensionless function 



C 2 (k, $) 



_ C 2 (KllSl ) 



Vfh 2r 



(A5) 



The resulting expression is then 

, , , d(d(> v n (Qv m ((>) 
c 2 (k, o) = I , X 



exp 



2ttV1 - k 2 

(( -C) 2 + 2(1 - K )(C + d)(C + S) 

2(1 -k 2 ) 



(A6) 



This is a general expression for arbitrary interaction po- 
tential; from now on we will make use of the specific form 
of v n (C)- Let us first consider the limiting cases k = 
and K = 1. First, for k = one gets 



c 2 (0,S) 



n 



k£ {m,n} 



— C fc e^ (c+l5)2/2 



(«5>i) T(n + l)r(m+ 1) _ 5 2 
~ 2?r S m+n + 2 6 

In the opposite limit «;—>!, 



p"(C+«5) 2 /2 



(<s»i) T(m + n + 1) _ 5 2 /2 



2tt 5-m+n+i 



(A7) 
(A8) 

(A9) 
(A10) 



Comparing the 5-dependence in eqs. (A?) and (A9) with 
the dependence of the mean distance d on the normal 
load, Eq. (||), - remember that we consider d to be a 
function of N - one finds to leading order 



c 2 (k,5) 



for k w 1 
for k « 0, 



(AH) 



where S is to be considered a function of Af/L 



D 



To make progress with Eq. (A6) for arbitrary k, one 
can introduce polar coordinates C = T cos (j), Q' = r sin <fi. 
The r-integration can be carried out exactly, leaving the 
one dimensional angular integral 

c 2 (k, 5) = T(p + 1) (1 - k 2 Y' 2 c-tt^ x (A12) 
"""" ^e^D_ {v+1){vm , 



cos sin 



where 



27r (1 - Ksin2</)) i V 



2 1 - n l + sm2<£ 2 

y (<p) = t— i r-^-r s , 

1 + k 1 — k sin 2<p 



(A13) 



and D_(„ + i)(y) with v = m + n + 1 is a parabolic cylin- 
der function. Equation ( |Al2j ) can easily be calculated 
numerically. For our present purpose it suffices however 
to write 

C2 ( K , S) = r( ^ +?1+1) f 2 (n) e-WW, (A14) 



where f 2 (n) only weakly depends on k ('weakly' com- 
pared to the leading exponential dependence). It is a 
monotonically growing function with 



h(K) 



r(»+i)r(m+i) i fnrK ._ n 

T{m+n+l) y^S iU1 K U (A15) 

1 for K = 1. 



8 



For large enough 6, we can neglect the dependence on 
/ 2 (k): The ratio e^ s ^ 2 /f 2 (0) is « 1/50 for (5 = 4, 
w 3 x 1(T 4 for 8 = 5 etc., with n = m = 2. 

Since we will primarily need the Fourier transform of 
the potential correlation function, we define 

<7 2 (p, p') = J d D x d D x' e ipx+ip ' x ' C 2 (x, x') (A16) 

and analyze its behaviour for large 5 such that e - * 2 / 2 <C 
1. Under this constraint, it is a good approximation to 
consider only the leading 6- and /{-dependence of c 2 (k, 6) , 
i.e. C2(k,S) ~ e -<5 "/( 1+K ', making an error of 0(1). For 
sake of simplicity, let us assume henceforth that both 
surface profiles obey the same distribution Kh{x) with 



I 2 
'o 



1/2 for|x|<cr 
for Ixl > a. 



Let 



~c 2 {p,p') = C 2 {p,p')/V 2 h 



2 7 2n 



(A17) 



(A18) 



Because C 2 (x, x') decays to a nonzero constant for |x| ^> 
<r, it is convenient to write c 2 (p, p') as the sum of four 
terms: 



c 2 5 P 5 P > + cf (p) 5p> + cf (p')Sp + c^(p, p' 



(A19) 



only the last two of which will contr ibu te in the integrals 
C 2 appears in (cf. for example Eq. (Bl)). The other two 
terms will be eliminated due to the 5 p -term being mul- 
tiplied with a power of p\ . 

Now, first consider Gaussian correlated surface profiles 
k(x) = e~ x I 2 (cf. the definition in Eq. ([!])). Perform- 
ing a saddle-point expansion about small x and x' in 
C 2 (x, x') and then Fourier transforming, cf(p) is found 
to be proportional to (a / 6) D e~ 2S / 3 for |p| < 6/a and 
exponentially damped for larger p. Correspondingly, 
c^(p,p') ~ (a/6) 2D e~ s ^ 2 for |p|, |p'| < 5/a, with expo- 
nential suppression for |p| or |p'| > 6/a. 

If the surface correlation function is taken to have a 
cusp in the xi-direction, the only - and important - 
difference arises in the large pi-behaviour of cf and c^. 
With the modification that 



(A20) 



leaving the dependence on the perpendicular spatial co- 
ordinates unmodified, we find 



<? 2 (px)~a h /{l+pial) 



(A21) 



with ah = 9a /26 2 , and a similar result for cl(pi,p'i). In 
the perpendicular direction p^, c 2 and c 2 are again ex- 
ponentially suppressed. 



APPENDIX B: DYNAMIC FRICTION FORCE 



We want to analyze the expression (|29j ) that gives the 
kinetic friction force density at velocity v in first non- 
vanishing order perturbation theory. Two generic scenar- 
ios will be considered, assuming in both cases that the 
two surface profiles obey identical distributions. First, 
we consider the case of short range correlated profiles, 
k{x) = e~ x / 2 , and second, profiles that are self-affine in 
the direction of F. Let us rewrite Eq. (p9|) as 



/ = 



G Q (p,qiv) iql C 2 (q,p - q), 



(Bl) 



p Jq 



with C 2 (p, p') as defined in Eq. (A16). Only the imagi- 
nary part of the propagator 



S(Go(p, <?!«)) 



7 2| p |2a + ( miV )2 



(B2) 



gives a contribution to the integral, the integral over the 
real part vanishes by symmetry. Now, C 2 (p,p') is of 
order 



V 2 h 2n e~ 52 / 2 ~ p fr for 1/L < |pO | < 1/a, 

(B3) 

with 6 as defined in App. |A| and from Eq. (p6|), and 
vanishes rapidly for |p^| > 1/a. This is the contribu- 
tion from c^(p,p'), cf. Eq. ( |A19 ). The second term in 



the denominator in (B2) serves as an infrared cutoff for 
the p-integration at p w a~ 1 (v/vo) 1 / a , giving rise to a 
velocity dependence 



<f>(v/vo) 



D-2a \Vq 

^ln(v/v ) 



D/a-1 

v ^ for D < 2a 



(B4) 



for D = 2a, 



hence confirming the discussion after Eq. (g^) for veloci- 
ties v <C vq. For t; > j)o, the 1/v- behaviour follows from 
the dom inance of the second term in the denominator of 
Eq. ( |B2[) , regardless of the dimension D and the value of 
a. The term cf(q)<5( p _ q ) gi ves an independent contribu- 
tion to the integral in Eq. (Bl) proportional to 



L 



7P' 



\r]vpi 



5f(p)- 



(B5) 



For small v <C vq, it is linear in v, and, apart from the 
velocity dependence, it is smaller by a factor of e~ s / 6 
compared to the leading term, so we can neglect it here. 

Now consider surface profiles that are ch aract erized by 
a correlation function modified via Eq. (A20), having 
a cusp in the direction of the applied force, and being 
analytic in directions perpendicular to it. In contrast 
to the previous situation, C 2 (q, p — q) is exponentially 
damped now only in the direction perpendicular to F, 
while the large p-behaviour of c 2 in the sliding direction 







(Pi 



3iJ 



2 q 1 2 , and 



p 1 2 for <§,. 



In the contribu- 



tion to Eq. (Bl) containing the term cj, performing the 
<7i-integration leads to 

2|Pl ° + ^ 2 C 2 (q^,p^-q^), (B6) 



2 7 cr 4 



(\p\ a +v) 2 



where v = v/vo<7 a . The remaining integrals can now eas- 
ily be carried out, yielding a result which, provided that 
D > a, is independent of v for v <C vq. It is of order 



1 



V 2 h 2n e- & I 2 



(D-a)~f<j 3 - a (D — a) L D 



Ffr- 



(B7) 



For a — 1, this result is multiplied with a factor of 
ln((r/a), where 1/a is the UV-cutoff of the pi-integra- 
tion, a setting the smallest length scale down to which 
the self-affinity of the surface profile holds, cf. Eq. (34). 
The result for the friction force in this case is indepen- 
dent of v in the regime (a/cr) 1 /"^ < t) C wo. For even 
lower velocities, the fact that the correlator in Eq. is 
analytic on length scales smaller than a comes into play 
again, leading to the same u-dependence of 4>(v/vo) in 
this regime as given by Eq. (B4). 

Again, there is an independent contribution from 
C2(q)<5( p _ q ) of the form Eq. (fB5|). This term has a de- 
pendence on the mean distance ~ e~ 2S / 3 and is hence 
proportional to ((a/L) D Af) 4 / 3 . For a = 2, this contribu- 
tion is negligible, but for a = 1, it suffers from a linear 
divergence in the UV-cutoff ~ a /a due to the large p\- 
behaviour of ■ We have checked that in higher orders 
perturbation theory no divergences in higher order of a /a 
appear, so this divergence will be compensated by the 5- 
dependence of this term for moderately large S. 
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